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Abstract 

We determine the neutrino mixing and mass parameters that are allowed by the 
Super-Kamiokande atmospheric neutrino data in a three-neutrino model with one 
mass-squared difference contributing to the oscillations. We find that although Ur 
oscillations are favored, i/^ — > Vf, oscillations with amplitude as large as 0.18 are allowed 
even after accounting for the limit from the CHOOZ reactor experiment. The range 
of allowed parameters permit observable ^ Vf. and i/g — oscillations in future 
long-baseline experiments. 



Introduction . It was suggested long ago [|l| that the atmospheric neutrino anomaly 
could be explained by the oscillation of muon neutrinos and antineutrinos into another 
neutrino species. This interpretation has been confirmed by the zenith angle dependence 
measured by the Super-Kamiokande (SuperK) experiment |0. Neutrino oscillations can also 
be invoked to separately explain the solar neutrino deficit 0, Q and the results of the LSND 
experiment 0. Because confirmation of the LSND results awaits future experiments and 
recent measurements in the KARMEN detector exclude part of the LSND allowed region [0, 
a conservative approach is to assume that oscillations need only account for the solar and 
atmospheric data. Then the two mass-squared difference scales in a three-neutrino model 
are sufficient to describe the data. Interest in the implications for models of the atmospheric 
neutrino anomaly has recently intensified 1^ . An attractive possibility is that both the 



atmospheric z/^ and solar Ue oscillate maximally or near-maximally at the Sm'^^^ and 5m^, 



scales, respectively 0, [10 



In this letter we use the recent Super-Kamiokande atmospheric neutrino data to deter- 
mine the allowed values for the general three-neutrino mixing matrix under the assumption 
that one mass-squared difference, 5m^j^, explains the atmospheric neutrino oscillations. In 
our scenario the other mass-squared difference, 5m <^ 6ml^^, can explain the solar neu- 
trino oscillations via either an MSW |rT| or vacuum long- wavelength scenario \T2\. We find 



that although pure Vr oscillations of atmospheric neutrinos are favored, there exist 

three-neutrino solutions with non-negligible z/^ oscillations, even after applying the 



constraints from the CHOOZ reactor experiment |T3[. Consequently z/^ <-> Ve and Ve 
oscillations may be observable in future long-baseline experiments. 

Oscillation probabilities. We begin our analysis with the survival probability for a given 



neutrino fiavor in a vacuum [|T4 



P(z/a ^ Z/a) = 1 - 4 ^ PajPak Sin^ Ajk , 



where 



Paj — ; (2) 

U is the neutrino mixing matrix (in the basis where the charged-lepton mass matrix is 
diagonal), A^-fc = 5m],^L/AE = 1.27((5m2^/eV^)(L/km)/(E/GeV), 6m]^. = m'j-ml, and the 
sum is over all j and k, subject to k < j. The matrix elements Uaj are the mixings between 
the fiavor {a = e, /i, r) and the mass (j = 1, 2, 3) eigenstates, and we assume without loss of 
generality that nii < m2 < m^. The solar oscillations are driven by IA21I = As„„ and the 
atmospheric oscillations are driven by IA31I ~ IA32I = Aatm ^ A^ 



The off-diagonal vacuum oscillation probabilities of this three-neutrino model are 
P{i^e ~^i^f.) = 4 PesP^s sin' Aatm ' 4.Re{U eiU^^U^^U ^2} sin' A,„„ - 2 J sin 2A,,„ , (3) 

P(z/e ~^Vr) = 4 Pe3^r3 siu' A^tm - 4Re{UelU:2U;^Ur2} slu' A,„„ + 2 J siu 2A,,„ , (4) 
Pil^^^^Ur) = AP^sPr3Sm'^Aatrn-4:Re{U^iU*2UrlUr2}sm^Asun-2Jsm2A,un, (5) 

where the CP-violating "Jarlskog invariant" is J = Y.k,'y^ijk^ai3'ylTT^{UaiUajUpJJpj} for 
any a, /3, z, and j (e.g., J = Im{Ue2U*^U*2U^3} for a = e, /5 = /i, i = 2, and j = 3). The CP- 
odd term changes sign under reversal of the oscillating flavors. We note that the CP-violating 
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probability at the atmospheric scale is suppressed to order Sml^^/ Sml^^, the leading term 
having cancelled in the sum over the two light-mass states; thus, Pii'a ^/s) = P{^f3 ~^ ^a) 
at the atmospheric scale. 

Fit to atmospheric neutrino data. For the L/E values of the atmospheric and long- 
baseline experiments, ^sun can be neglected, and the vacuum oscillation probabilities become 
simply 



PU 



1-4P^3(1 



P,,^) sin^ A„ 



P(iy, 



^3 ; L^atm 1 

l-4Pe3(l-Pe3)sin2A„ 
4Pe3P/.3 sin^A 

4P^3(1-Pe3 



^atm J 



4P.,fl-P.,-P^3)sin2A,j„, 



P^3) Sin^ Aatm 



(6) 
(7) 
(8) 
(9) 
(10) 

In Eqs. (|]) and ( p!OD we have used the unitarity condition P^^ = 1 — Pgs — P^3. Thus for 
oscillations at the atmospheric scale there are only two independent mixing matrix parame- 
ters, e.g., Pe3 and P^3, that are relevant. All predictions for atmospheric and long-baseline 
experiments are completely determined by the three parameters 5m^^^, Pe3, and P^3. We 
define the oscillation amplitudes A'^f^^, v4[jf^, Aatm: A^^, and A^lm^ as the coefficients of 
the sin^ Aatm terms in Eqs. (|6|)- ([ToD , respectively. The parameters Pe3 and P^3 can then be 
determined from the atmospheric neutrino data by the relations 



a 



{S) Aatm) + {S) Aatm 



(11) 



and 



a 



'1 



(S) A"'^ 



atm, 



+ r ^ {S) Aatm 



(12) 



where N" and A*"" are the expected numbers of atmospheric e and fi events, respectively, 
r = N°/N°, {S) is sin^ Aatm appropriately averaged, and a is the overall neutrino flux 
normalization, which we allow to vary following the SuperK analysis 0. SuperK reports 
N^/N° and N^/N" from a 535 day exposure for eight different L/E bins [Q. The data were 
obtained by inferring an L/E value for each event from the zenith angle di and energy of the 
observed charged lepton Ei and comparing it to expectations from a monte carlo simulation 
based on the atmospheric neutrino spectrum 



16 folded with the differential cross section. 



Due to the fact that the charged lepton energy and direction in general differ from 
the corresponding values for the incident neutrino (or antineutrino), the L/E distribution 
involves substantial smearing. We estimate this smearing by a monte carlo integration over 
the neutrino angle and energy spectrum |]l^ weighted by the differential cross section. We 
generate events with E,y and 6^, and determine the corresponding Eg and 6e for the charged 
lepton. We bin the events in L/E,^, using to determine L and an estimated neutrino energy 
inferred from the average ratio of lepton momentum to neutrino energy, E^^^ = Eg [E^/ Ei) , 
analogous to the SuperK analysis We then calculate a value for (sin^ Aatm) for each 



L/E bin for a given value of (5m^^^. Finally we make a fit to Eqs. (|I^) and (|T2p to determine 



Pe3^ P 



m3: 



^^Itm^ and a 



18 1 . Without loss of generality we take SrUat^ to be positive. 



Our best fit values for the four parameters are 



6m] 



atm 



2.8 X 10"^ eV 



(13) 
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Pes 
P, 



0.00, 

;.3 = 0.50, 

a = 1.16, 



(14) 
(15) 
(16) 



with Xmin = 7.1 for 12 degrees of freedom. This best fit is close to the result of the SuperK 
simulation that assumed only two-neutrino oscillations. In Fig. la we show the allowed 
region for P^3 versus Pgs for a = 1.16 and Sm^^j^ = 2.8 x 10^'^ eV^. Although Pes = 
is favored, small nonzero values are allowed, which permit some ^ z/g and z/g — > Ut 
oscillations of atmospheric neutrinos. In Fig. lb we show the allowed region for the overall 
flux normalization a versus Sml^^ for P^s = and P^s = 0.50. 



Another limit on Pes comes from the CHOOZ reactor experiment [|T^] that measures z/g 
disappearance 



A 



atm 



4Pe3(l-Pg,) <0.2 



eSj 



;i7) 



which applies for Sml^^ ^ 2 x 10 ^ eV^. 



varies with 5m^^„, and 



for SiTT'ltm 



eV and a = 1.16, 



^ ^ xw V. V . The exact limit on Aatm 
< 10-3 eV^ there is no limit at all. For 5m^j^ = 2.8 x 10 ^, ^^^^ ^ — x.^^, 
Pe3 is constrained to be less than 0.04. The result of this additional constraint is shown in 
Fig. la. In Fig. Ic we show the effect of the CHOOZ constraint on the allowed region of Pes 
versus Sm^^^ for a = 1.16 and P^^a = 0.50. 

Varying over the entire parameter space, and imposing the CHOOZ constraint, the ranges 
of allowed values at 68% (95%) C.L. are 



0.8(0.5) < 6ml,J{10~'eV')< 7.9 (10.0) 

0.00(0.00) < Pe3 < 0.05 (0.08) , 

0.29(0.25) < P^3 < 0.71 (0.75) , 

1.07(1.04) < a < 1.24(1.28) . 



(18) 
(19) 
(20) 
(21) 



The allowed ranges of some related oscillation parameters are 



0.29 (0.25) 
0.83 (0.75) 
0.00 (0.00) 
0.00 (0.00) 
0.00 (0.00) 
0.82 (0.74) 



< 
< 
< 
< 
< 
< 



Pr3 < 0.71 (0.75) 
A^i< 1.00 (1.00) 
A:L< 0.18 (0.29) 
A::L< 0.11 (0.18) 
A:I^< 0.10 (0.16) 
A:^L< 1.00 (1.00) 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



fir oscillations are strongly favored 



are allowed with amplitude as 



Although 
large as 0.18. 

We do not consider matter effects in the above analysis. For the Srn^^j^ favored by our 
fit, matter effects are small for the sub-GeV neutrinos that constitute most of the data. 
Also, as evidenced by our fit, the dominant oscillation is z/^ — *• u^, which is not affected by 
matter. However, matter effects could be important for neutrinos with smaller 6m'^^^/E, i.e., 
especially for multi-GeV data in solutions with Sml^j^ ^ 10^^ eV^, which could be observable 
as more data becomes available at higher energies. 
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New physics predictions. The new physics predicted when Pgs 7^ is z/e ^ z/,- oscilla- 
tions with leading probability (i.e., the /S.atm terms) given by Eq. (H). From Eq. (|^) the 
maximal amplitude A^J^^ for these oscillations is 0.16 at 95% C.L. Figure 2 shows the al- 
lowed values of A^J^^ versus v4^^^ for 5rn? = 2.8 x 10~^ eV^ and a = 1.16; the effect of 
the CHOOZ constraint is also shown. The z/g Ur oscillations could be observed by long- 
baseline neutrino experiments with proposed high intensity muon sources [|l^, 20, which 
can also make precise measurements of z/^ ^ z/g and z/^ —>■ Vj. oscillations. Sensitivity to 



'^'atmi.^^^tm/ Y > 2.5 X 10 ^ is cxpcctcd for the parameter ranges of interest here; 



for 5m? = 2.8 x 10"^ eV^, A'^^^ could be measured down to 3 x 10~^. The measurement 
of z/e Vr and ^ Vf, oscillations in such a long-baseline experiment would test the 
three-neutrino model. 

Discussion . Further measurements of atmospheric neutrinos will more precisely deter- 
mine the parameters Pes, P^s, and 5ml^„^. The MINOS [||], K2K [||], and ICARUS |4 



experiments are expected to test for — > z/g and z/^ — > Ur oscillations for 6ml^^ > 10~^ eV^. 
Together these measurements could put strong limits on Pgs, which governs the z/g Vr 
oscillations that could be seen in future long-baseline experiments such as those utilizing 



a muon storage ring at Fermilab |T9[. Full three-neutrino fits including the solar neutrino 
data 1^ can then determine one of the remaining two independent parameters in the mixing 
matrix, e.g., Pei, using the z/g survival probability for Aatm 3> 1 

P(z/g ^ i/g) = 1 - 2Pg3(l - Pes) - 4Pei(l - Pgi - Pgs) siu^ A,„, . (28) 

The considerations in this paper can also be extended to a four-neutrino model ||21|| . 
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Figure 1: Allowed regions of the four-parameter space ((5m^ 
line) and 95% C.L. (solid line) for (a) P^s versus Pes with 5m^j„ 



a 



x-es, -r;,3, a) at 68% (dashed 
= 2.8 X 10-3 eV^ and 
1.16, (b) versus a with Pe3 = and P^3 = 0.50, and (c) 5m^^„ versus Pes with 



-P/iS = 0.50 and a = 1.16. The best fit values are indicated by the filled circles. The CHOOZ 
constraint |13[, shown by a dotted line, excludes the shaded region. 
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0.2 




Figure 2: Allowed regions at 68% (dashed line) and 95% CL. (solid line) for the oscillation 
amplitudes Agr versus with 6m'^ = 2.8 x 10^'^ eV^ and a = 1.16. The best fit values are 



indicated by the filled circle. The CHOOZ constraint 
the shaded region. 



shown by a dotted line, excludes 
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